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Motivation 


□  Microstructure  and  properties: 

>  Material  properties  are  sensitive  to 
microscale  structure. 

>  Microstructure  are  inherently 
random  but  correlated. 

>  Microstructure  variation  induces 
property  variability. 

□  Importance  of  uncertainty  analysis: 

>  Assess  product  and  process 
reliability. 

>  Estimate  confidence  level  in  model 
predictions. 

>  Identify  key  sources  of  randomness. 

>  Provide  robust  design  solutions. 

□  Focus:  Polycrystalline  materials 

□  Solution  strategy: 

>  Construct  stochastic  input  model. 

>  Solve  the  stochastic  partial 
differential  equations  (SPDEs). 

Cornell  University 

Materials  Process  Design  and  Control  Laboratory 


3 


Convex  hull 


/-*  .  Equivalent  strain  , 

Stress-strain 


Property 

200 

=150 


distribution 


Primary  Development 


□  Employ  model  reduction  techniques  to  construct  the  reduced-order 
surrogate  model  of  random  microstructures  based  on  given  samples. 

□  Develop  physics-based  deterministic  solvers  to  estimate  mechanical 
properties/responses  of  polycrystalline  materials  based  on  the 
interrogation  of  microstructures. 

□  Investigate  property/response  variability  of  polycrystalline  materials 
using  efficient  stochastic  simulation  methods. 
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Problem  Description 


>  Problem  definition 

Given: 

■  Grain  size  snapshots  constrained  by  moments  (mean  size,  standard  deviation, 
higher-order  moments). 

■  Texture  snapshots  from  random  process. 

Goal: 

■  The  variability  in  material  properties  and  responses 

>  Methodologies 

■  Model  reduction  to  reduce  the  complexity 
of  stochastic  input 

S  Nonlinear  Model  Reduction  (isomap)  to 
reduce  grain  size  space 
S  Karhunen-Loeve  Expansion  to  reduce 
texture  space 

■  Crystal  plasticity  Taylor  model  for  property 
estimation. 

■  Adaptive  sparse  grid  collocation  to  solve 
stochastic  partial  differential  equations 

(Z.  Li,  B.  Wen  and  N.  Zabaras,  2010) 
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Microstructure  Representation 


Microstructures  obtained  from  certain  process  satisfy  some  specific  experimentally 
determined  statistics  of  grain  size  distribution. 

Each  microstructure  that  satisfies  the  given 
statistics  of  the  grain  size  distribution  is  a 
point  that  lies  on  a  manifold  embedded  in  a 
high-dimensional  space. 


For  microstructures  having  the  same  mean  grain 
size,  a  “sorted  grain  size  vector”,  whose 
dimension  is  invariant,  can  be  used  to  carry  the 
grain  size  information. 


-15  15 


The  difference  between  two 
microstructures  is  conveniently 
measured  by  Euclidean  distance. 


D(A,B)  = 


Microstructure 

construction 


'±(GS?-GS?)2 


,1/2 


V  «•= 1 
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Nonlinear  Model  Reduction  on  Grain  Size  Feature 


Given  a  set  of  N  unordered  points  belonging  to  a  manifold  M  embedded  in  a  high¬ 
dimensional  space  M1,  find  a  low-dimensional  region  A  c  HP  that  parameterizes 
M,  where  d  «  n. 

1 )  Geometry  can  be  preserved  if  the  distances  between  the  points  are  preserved  - 
Isometric  mapping. 

2)  The  geometry  of  the  manifold  is  reflected  in  the  geodesic  distance  between  point. 


Algorithm: 

1)  Compute  the  low-dimensional  representation  of  a  set  of  N  unordered  sample  points 
belonging  to  a  high-dimensional  space 


Given  N 

Compute 

Perform  MDS 

N  points  in  a  low 

unordered 

pairwise 

on  this 

dimensional 

samples 

geodesic 

distance 

distance 

matrix 

space 

2)  For  an  arbitrary  point  y  e  A  must  find  the 

corresponding  point  x  e  M.  Compute  the  A  C  HP  Mr  HP 

mapping  from  A-^M  based  on  k-nearest  J  c_ 

neighbors.  . ***. 

— —  * 

(B.  Ganapathysubramanian  and  N.  Zabaras,  2008) 

foJiiOl  Cornell  University  8 

Materials  Process  Design  and  Control  Laboratory 


Texture  Modeling 


The  properties  of  a  polycrystalline  microstructure  are  highly  dependent  on  its 


,  .  h 

,  *  ^=> 


b 


texture:  orientation  distribution  of  grains. 

Orientation  representation:  Rodrigues  parameters 

6  6  6 

rx  =  wx  tan  — ,  r2  =  w2  tan  — ,  r3  =  w3  tan  — 

2  2  2 

Texture  representation:  Grain  orientation  vector,  Orientation  Distribution 

Function  (ODF) 

A  discrete  form :  x(r)  =  jr/ ,r2  ,r31,...,rln /2 , r3 } 


12  slip  systems  in  FCC 


A n  w 


ODF  in  RF  fundamental 
zone  of  FCC  crystal 


n*w:(llO) 

Kca, :  (1 1 1} 


Orientation  dependence  of  slip  system  (anisotropy  in  crystalline  materials) 


=  R'mL/ 


where  Ry  =- — t — _ |l(l  —  r1  -rJ^  +  2(rJ  ®rJ  —  Ixr'^j 
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Generation  of  Initial  Texture  Samples 


Initial  texture  samples  can  be  obtained  by  a  sequence  of  random  processing 
simulations  with  various  deformation  rate 
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co1,co2,...,cds  are  random  coefficients  corresponding  to  tension/compression, 
plain  strain  compression,  shear  and  rotation. 

The  slip  systems  are  updated  during  deformation  as 

=Fe(0i< 
n  *  =F  e~T{t)K 

Therefore  the  new  orientations  can  be  recovered  from  the  rotation  part  of  the 
elastic  deformation  gradient. 
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Karhunen-Loeve  Expansion  on  Texture  Samples 


Deterministic  texture 


Initial  texture  samples 


Random  process 
controlled  by 

|  CQ^ ,  CO  2  5  •  •  •  5  (Og  | 


•  •  • 


Given  N  texture  examples,  construct  covariance  matrix  of  these  samples 


t  = 


The  truncated  Karhunen-Loeve  Expansion  of  a  random  vector  t  is 

d 

T(r,to)  =  f(r,co  )  + 

i— 1 

where  are  the  /th  eigenvector  and  eigenvalue  of  C  ,  respectively.  {^O)}  are 
a  set  of  uncorrelated  random  variables  satisfying 

(«*)) = E{n, ,(®)7y(®))=4> 

Texture  random  field  thus  transformed  to  low-dimensional  space  r/eRd 
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Maximum  Entropy  Estimation  of  the  Distribution  of  r/ 


To  sample  new  texture,  we  can  sample  77  instead,  and  then  transform  it  back  to  the 
texture  space.  The  distribution  of  77  is  needed. 

Maximum  Entropy  Estimation  (MaxEnt):  amongst  the  probability  distributions  that 
satisfy  our  incomplete  information  about  the  system,  the  probability  distribution  that 
maximizes  entropy  is  the  least-biased  estimate  that  can  be  made.  In  agreement 
with  everything  that  is  known  but  carefully  avoids  anything  that  is  unknown. 

The  form  of  MaxEnt  distribution  is 

~YuAnf(x)  " 

p  «= 1  r  ~YjXnfix) 

p*(x)  = - ,  Z-  J  e  n=1  dx 

Z 

which  maximize  the  entropy 

M 

H(p)  =  ~2  P(x,  )  log  (p{xi )) 

and  satisfies  constraints 

ev.W)=m„  n  =  1,2,... 

When  the  uncorrelated  constraints  are  satisfied,  the  MaxEnt  distribution  is  a 
standard  Gaussian  distribution 

//-AW) 
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Inverse  CDF  Transformation 


Uncorrelated  Gaussian  random  variables  are  independent  in  current  case.  To 
employ  Sparse  Grid  Collocation  method,  Gaussian  distribution  needs  to  be 
transformed  to  a  uniform  hypercube  [0,1  ]d. 


The  cumulative  distribution  function  (CDF)  for  standard  Gaussian  is 

l 


%M)= 


l  +  erf 


7,  ^ 


V 


V2 


which  is  uniformly  distributed  in  [0,1]. 

Given  a  point  in  the  hypercube  ?e[0, i]d,  we  can  also  find  a  corresponding  point 
in  the  original  distribution  by 

77.  = 


This  process  transforms  a  node  in  sparse  grid  back  to  a  point  in  Gaussian 
distribution,  and  it  can  be  further  recovered  to  a  texture  realization 
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Adaptive  Sparse  Grid  Collocation 


Sparse  grid  collocation  is  an  effective  method  to  solve  SPDEs.  It  approximates  the 
multi-dimensional  stochastic  space  using  interpolating  functions  on  a  set  of 
collocation  points.  The  collocation  method  collapses  the  multi-dimensional  problem 
to  solving  M  (M  is  the  number  of  collocation  points)  deterministic  problems. 


The  interested  function  can  be  approximated  by 


u(x,%(G>))  =  Yj  Z  0)\{x)a\i^{co)) 

I _ I  j  I _ I  I _ I 

4  4  4 

Stochastic  Hierarchical  Multi-linear 
process  surplus  basis  functions 

The  mean  of  the  random  solution  is  evaluated  as 


E(u(t))  =  YdYJCo\(x)-\  a\&)d% 

|i|*9  J' 


In  the  context  of  adaptivity,  new  support 
nodes  are  added  to  the  hypercube  only  if  the 
error  indicator  is  larger  than  a  threshold  z\ 


a>\(x)-\La\(%)dt, 


\En. 


i  n-d  - 1 


>  £ 
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Deterministic  Solver 


Material:  FCC  nickel 

Deterministic  solver:  rate-independent  crystal  plasticity  with  Taylor  homogenization. 
Decomposition  of  deformation  gradient:  F  =  F^F^ 

Hardening  law:  k-k^  =  ajubjp 

where  dislocation  density  changing  rate  =  2 

K 

If  the  resolved  shear  stress  is  larger  than  the  slip  resistance  ra  >  Ka,  the  slip 
system  is  active. 


l-  +  kiy[p-k2p 


L§b 


r 


Upon  calculating  the  incremental  shear  strain,  elastic  and  plastic  deformation 
gradient  can  be  updated  and  Cauchy  stress  is  computed 

T  =  Fe  { [det  Fe  ]_1  T*  }¥eT  with  T*=CeEe 


Homogenized  effective  stress  and  strain 

T  =  (T)  =  — f  TdV 

\  /  V  JV 


3-  - 

T'«T' 

2 


where 


£*=lhD-mt 


V 

"  ’  D  =  (D)  =  — f  DdF 

\  /  y  Jv 


(L.  Anand  and  M.  Kothari,  1996;  A.  Acharya  and  A.J.  Beaudoin,  2000) 
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Investigating  Property  Variability  Due  to  Microstructure  Uncertainties 


Property  variability 


Property  convex  hull 


Property  distribution 


Microstructure  space 
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Reconstruct  microstructures 


Low-dimensional  space 


START:  Extract  microstructural  features 
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Numerical  Examples 


Microstructure  domain:  1mm  x  1mm  x  1mm  cube  containing  54  grains. 


Examine  the  effective  stress  variability  of  polycrystalline 
microstructures  satisfying  different  constraints  subjected  to 
compression. 


L  =  0.002  sec  1 


0.5  0  0 

0  0.5  0 
0  0-1 


Grain  size  sample  constraints: 

(1)  Mean  volume  0.0185  mm3 

(2)  2nd  moment  3.704x1 0_4mm6 

(3)  3rd  moment  8.637x1 0_6mm9 


Texture  generation: 
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with  g\,cd2  ~  [-0.002,  0.002] sec  1 

Run  the  simulation  for  500  seconds. 
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Numerical  Examples 


JQd 


Elastic  Properties  Variability 


Shear  modulus  distribution 
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Overview 


□  Nickel-based  superalloys  are  widely  used  in  components  working  in 
harsh  environment,  e.g.  disks  and  blades  of  turbo  engines. 

□  They  have  high  resistance  to  creep  and  fatigue  at  high  temperature. 

□  The  performance  of  nickel-based  superalloys  depends  on  their 
underlying  microstructures:  grain  topology,  texture,  and  volume 
fraction  of  y’  phase. 

□  We  are  interested  in  the  variability  of  fatigue  properties  of  superalloy 
microstructures  due  to  uncertainties  in  texture,  grain  size,  and 
volume  fraction  of  y\ 

□  Principal  component  analysis  (PCA)  based  model  reduction 
techniques  are  adopted  for  reducing  the  complexity  of  the  input 
space. 
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Problem  Definition 


>  Problem  definition 


Given: 

■  Grain  size  snapshots  constrained  by  the  mean 
size  and  sampled  from  lognormal  distribution. 

■  Texture  snapshots  from  random  process. 

■  Volume  fraction  of  y’  precipitates. 

Goal: 

■  The  variability  of  fatigue  resistance  due  to 
microstructure  uncertainties. 

■  Find  the  feature(s)  dominate  the  variability  of 
fatigue  resistance 

>  Methodologies 

■  PCA  and  KPCA  model  reduction  to 
reduce  the  complexity  of  stochastic  input. 

■  PCE  to  map  reduced  coordinates  to  a 
known  distribution. 

■  Adaptive  sparse  grid  and  Monte  Carlo 
collocation  to  solve  stochastic  partial 
differential  equations 


maxPFS 


(Bin  Wen  and  N.  Zabaras,  2012) 
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_ Stochastic  Input  Space _ 

□  A  probability  space  (£2,F,P)  with  sample  space  Q  ,  corresponding  to 
all  microstructures  resulted  from  certain  random  process,  Fc2°  , 
being  the  o-algebra  of  subsets  inf2  ,  and [0,1]  ,  as  the 
probability  measure. 

□  Each  sample  co  e  Q  is  a  random  microstructure  which  can  be 
described  by  a  discretized  representation  y(<x>)  =  (yl,...,yM)T  YeRM 

□  The  input  space  information  is  only  given  as  a  set  of  samples 

y,.,z'  =  l,...,2V 

□  The  variability  of  microstructure  sensitive  property  A  =  A( y)  of  the 
given  input  space  is  interested. 

□  It  is  necessary  to  construct  the  input  space  that  has  the  same 
statistical  properties  with  the  given  samples.  Since  the  input  space  is 
high-dimensional,  model  reduction  is  adopted  to  construct  a  low¬ 
dimensional  surrogate  space. 

[HBB  Cornell  University  23 
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PCA/KPCA 


□  In  general,  the  components  of  stochastic  input  yeMM  are  not 
linearly  related. 

□  Direct  linear  PCA  attempting  to  fit  a  linear  surface  such  that 
reconstruction  error  is  minimized  is  not  very  appropriate  to  find  the 
surrogate  space. 

□  Solution:  nonlinearly  map  the  given  samples  to  a  feature  space  F, 
where  PCA  performs  well. 

y  x  =  <3>(y) 


Linear  PCA  Kernel  PCA 

5>:R"  ^F,y->x 


□  If  PCA,  ®(y)  =  y 

(B.  Schlkopf,  A.  Smola  and  K.-R.  Muller,  1998;  X.  Ma  and  N.  Zabaras,  2011) 
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PCA/KPCA  Formulation 


Covariance  matrix  in  the  feature  space: 


\  N  ~ 
^  i= 1 


O  =  O(y)  -O 


Eigenvalue  problem: 
KL  expansion: 


CV  =  XV 


Ms 


Me 


®(y) = Z  yfivn  + ® = 2>,  v,  +  ® 


/=! 


/=! 


Multiply  with  O  and  reform  the  eigenvalue  problem: 

KU  =  AU  where  K:^..  =  (4>(y,)-®(y,-)) 

U  [tt^  5  ...  5  tt  jy  ] 

A  =  diag^,...,AN) 

The  eigenbasis  of  the  covariance  matrix  c  is  projected  to  the  given 
samples  by  N 

v(=Z«#®(yy)=Z-7f®(yy) 

j=i  j=i  y/lz- 
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PCA/KPCA  Formulation 


TV 


TV 


Any  sample  in  F:  ®(y)  =  5>,.v,.  +  +  ®  =  XA®(y,-) 


i= 1 


i= 1 


z=l 
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where  p  =  AZ+ — 1  a  =  hu  = 

TV 


f  1  ^ 

I - llr 


V 


TV 


[«!,.. .,or]  Z  , . . .  ?  ] 


y 


Reduced  representation: 

Z,  =  v,  -®(y)  =  |;a(,(®(y)-4Kyy))  =  |;a#*(y,yy)  =  afk,  =  afHk,  --jUfHKl, 


Gaussian  Kernel:  kv  =  A(y,,y/)  =  exp 


f 

yi-yj 

2  \ 

V 

2  cr2 

J 

InPCA  £(y,,y,)  =  y, -y,- 


Only  the  dot  products  of  vectors  in  the  feature 
space  are  required,  while  the  explicit  calculation 
of  the  map  is  not  necessary. 


Construct  reduced-order  space  by  initial  samples 


Random  sample^  in  the  reduce  space  (constructed  byz,),  we  can  obtain 
its  corresponding  high-dimensional  counterpart. 
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Polynomial  Chaos  Expansion  (PCE) 


Map  the  reduced  space  to  a  known  distribution. 

/  \  ^  Each  independent  random  variable  £  can 

ih  ’  ’  be  expanded  on  to  an  one-dimensional 

polynomial  chaos  basis  Vj  of  degree  p. 

Gaussian-Hermite  PCs 


n 


E\4,vM ?<)]  i  r  e  ,  ^  -f  j 

=1=  rr  ,1  4iVj(fh)e  dlln  « =!,•••,'',  J=0,...,p 

■J2n  ir~° 


Wm)) 


KJ 


Uniform-Legendre  PCs 


r,v  =  — r 


E\Z,wM)\  2  j  +  \ 


[vM id) 


W%iY,(rh)& 1„  i  =  h—,r,  j  =  0,...,p 


To  compute  the  integral,  a  map  between  £  and  n,  is  needed. 

Mapping  through  CDF 

6  =ri07,-)>  r(SFfi-‘o^ 

where  and  T/;  are  the  CDFs  of  the  two  random  variables,  respectively. 

(D.  Xiu  and  G.E.  Karniadakis,  2002;  G.  Stefanou,  A.  Nouy  and  A.  Clement,  2009) 
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Pre-imaging 

□  Sample  11 ,  and  find  the  input  realization  y,  approximately. 


Draw  random  sample  if 

I 

Map  to  %  PCE 

I 

Construct  high-D 
representation  O(y) 

J 

Recover  microstructure  feature  y  Pre-imaging 

y  *  J_  Gauss  kernel: 

Z"-1  d.  y 


®(y) 


PCA/KPCA 


> 

/ 

k 

R2 

V 

s 

\  " 

\ 

■ 

cK 


* 


®  (Y)? 


Linear  PCA 


Kernel  PCA 


l  +  prKp-2prk 


y* 


K-nearest  neighborly  ~  y  = 


K  1 

^i=lJ 


d,2(y,yt)  =  ||y  -y,|f  =  -2  V1  ln(l  -0.542(®(y),®(y,.») 

“  -2ct2  ln(l  -  0.5rf,2(®(y),®(y())) 
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Crystal  Plasticity  Constitutive  Model 


□  The  superalloy  microstructure  is  modeled  as  a  homogenized  single 
crystal  with  effective  properties. 

Secondary 
Tertiary  y 


Explicit  structure 


■  £ 


Homogenized  model 


□  Taylor  model:  all  grains  have  the  same  deformation.  No  realistic 
microstructure  is  required.  Computationally  efficient. 

□  Finite  element  model:  the  boundary  nodes  have  the  same 
deformation  while  heterogeneous  deformation  is  allowed  within  the 
microstructure.  Realistic  microstructure.  More  accurate  but  less 
efficient. 


Cornell  University 

Materials  Process  Design  and  Control  Laboratory 


29 


Homogenized  Constitutive  Model 


□  The  model  is  for  polycrystalline  INI 00.  All  grains  are  assumed  to  be 
the  mixture  of  y  matrix  and  y’  precipitates.  The  same  set  of 
constitutive  equations  apply  to  all  grains. 

□  Three  types  of  y’  precipitates:  primary,  secondary,  and  tertiary. 

□  The  effect  of  y’  precipitates  is  taken  into  consider  as  constitutive 
parameters  (volume  fraction,  mean  particle  size,  etc.). 

□  Rate  dependent  flow  rule: 


■x 


O) 


D 


(a) 


+n 


(a) 


T(a)  _  %(a) 


D 


(a) 


□  Hardening  law:  ep  =  +  a.njbj 

where  a,  =  (o.  1  -  0.68 fpl  + 1 . 1  f'p\ ) 


Pmix  =  (fpl  +  fpl  +  fpl)Pf  +  frMr 

fp 3  J  oct  octahedral  slip  system 

fp3+fr  [cub  cube  slip  system 


K 


(a) 

0  ,oct 


K 


(a) 

0  .cub 


+  Ylb 


”k 


Vi 


"pi- 


fpi 

w — —  +c 


d, 


p  2 


w- 


f 


p  2 


+  c 


P3 


o  +  c  d 

3  gr  gr 


-0.5 


Y’  hardening 


Grain  size 
effect 
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APB 


APB  ref 


(C.P.  Przybyla  and  D.L.  McDowell,  2010) 
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Homogenized  Constitutive  Model 


□  Dislocation  density  evolution: 


Pi)=h 


k* 


bd 


Seff 


r 


(a) 


d 


f  2  V* 


Seff 


\d2 s  j 


maximum  range  of  cyclic 
plastic  shear  strain 


□  Back  stress  evolution: 

=CX[WrniXbif¥  ^(r{a)  -Z[a))-Z[a))j 

C  =  123.93  -  433.98 f2  +  384.06 f'l  d 


ArL  = 

u 


f 


maxfr*  )-min(^) 

eye  x  7  eye  x  7 


(«) 


r 


170.A 


Z0  +  kXx 


P[a) 


□  Fatigue  indicator  parameters  (FIPs)  are 

employed  to  measure  the  fatigue  properties 
of  nickel-based  superalloy.  r  . _ 

□  Cumulative  plastic  strain  per  cycle  ^=UI^=u  jDf:DPdt 

□  Cumulative  net  plastic  shear  strain  pr  =  max(j 

□  The  Fatemi-Socie  parameter  pfs=— 


\ 


(ht 


■rS 


S^nydt 


1 +k 


□  Maximum  range  of  cyclic  plastic  shear  strain  pmps  =— 

(M.  Shenoy,  J.  Zhang,  and  D.L.  McDowell,  2007) 


P 

max 
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Value 


□ 


□ 


Numerical  Example:  Taylor  Model  vs.  FEM 

Polycrystalline  microstructures  with  homogenized  grains 
implemented  in  the  Taylor  model  and  FEM. 

Initial  input  samples:  1000  randomly  generated  microstructures 
composed  of  54  grains  represented  by  grain  size  and  texture 
features. 


=0.42,/„,  =0.11 


P  2 


Taylor 


-0.70%  ~  0.70% 

0.35 

0.3 

*^Pcyc 

^^Pr 

0.25 

-£-PFS 

0.2 

)(  Pmps 

0.15 

0.1 

0.05 

0  I 

s  -  0.001s 


-1 


0.5 


i 


Strain  ~ 

' 


maxP  =1.51x10  ,maxP  =1.12x10  ,maxi%  =6.50x10  ,maxP  =5.98x10 
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Generation  of  Initial  Samples 


□  Initial  texture:  obtained  by  a  sequence  of  random  processing  simulations 
with  various  deformation  rate 


1 

0 

0 

0 

0 

0 

0 

-1 

0 

II 

£ 

0 

-0.5 

0 

+  co2 

0 

1 

0 

+  co3 

1 

0 

0 

0 

0 

-0.5 

0 

0 

-1 

0 

0 

0 

a>vco2,co3  are  random  coefficients  corresponding  to  tension/compression, 
rotation,  and  shear,  respectively. 

The  slip  systems  are  updated  during  deformation  as 

n,  =F  (0^0 


□  Initial  grain  sizes:  obtained  by  sampling  in  a  lognormal  distribution.  The 
mean  grain  size  is  fixed  at  0.0265  mm. 


ln(  <4  )  -  ) 


dgr  =  exp(//  +  a2  /  2)  =  0.0265 
cr  =  0.025 


□  The  above  process  is  only  for  generating  data  having  inherent 
correlation.  After  that,  the  only  accessible  information  of  the  input  is  the 
data.  No  knowledge  of  how  they  are  generated  is  assumed. 

Cornell  University  33 

Materials  Process  Design  and  Control  Laboratory 


PCA/KPCA  on  Texture 


4,  =  0,  /„2  =  0.42,  /„3  =  0. 1 1  E  =  -0.7%  ~  0.7% 


/*2 


>3 


s  =  0.00ls 


-1 


reduced  from  162-dim  to  4dim)  4  components  in  KPCA:81.5% 
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PC  expansion  on  reduced  random  variables 


£ =!>.#,  to)-  /=1> 


7=0 


r 


Uniform-Legend  re  PCs  give  better  representation  of  the  reduced  random  variables 
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FIPs 


MC-PCA 


max  Pcyc 


max  Pcyc 


MC-KPCA 


Cornell  University 

Materials  Process  Design  and  Control  Laboratory 


36 


FIPs 


ASGC-PCA 


ASGC-KPCA 


max  P-g 


ASGC:  level  8,  converge  to  err=0.0001 
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Convex  Hul 


PCA 


KPCA  s 


Some  of  the  tails  predicted  by  ASGC  are  not  captured  by  MC. 
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Grain  Size  Variation 


□  Fix  texture  and  volume  fractions  of  y’  particles,  while  vary  the  grain  sizes 
(sampled  from  the  lognormal  distribution). 

□  The  reduced  dimensionality  is  taken  to  be  10,  which  captures  91.2%  of  the 
total  energy  by  PCA,  and  89.2%  of  the  total  energy  by  KPCA. 

max  Pcyc  x  max  ft  x  iq*  max  PFS  x  10*  max  PmPs 


□  Distributions  of  maximum  FIPs  extracted  from  the  1000  initial  sample 
microstructures:  the  texture  dominate  the  FIPs  variability. 

max  Pcyc  x  10*  max  Pr  maxPFS  max  Pmps 


Combined  Volume  Fraction  and  Texture  Variation 


L 1  =0,/,2  ~[/(0.3, 0.5), fDi  ~ C/(0.11, 0.1 4) 


>2 


^3 


Texture  dim =4 


MC-PCA 


x  10"* 

r__ 

0  0125  0  013  0  0135  0  0U  0  0145  0  015  0  0155  0.016 

maxP-,. 


ASGC-PCA  ‘ 
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Outline 


□  Introduction  and  motivation 

□  Uncertainty  quantification  at  a  single  material  point 

>  Investigating  mechanical  response  variability  of  single-phase 
polycrystalline  microstructures 

>  Investigating  variability  of  fatigue  indicator  parameters  of  two- 
phase  nickel-based  superalloy  microstructures 

□  Uncertainty  quantification  of  multiscale  deformation  process 

□  An  efficient  image-based  method  for  modeling  the  elasto-viscoplastic 
behavior  of  realistic  polycrystalline  microstructures 

□  Conclusion  and  future  research 
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Multiscale  Model  Reduction 


□  Motivation:  Microstructure  features  are  location  specific.  Traditional 
construction  of  reduced-order  stochastic  models  for  one  point  cannot 
see  the  microstructure  correlation  between  points  in  the  macroscale. 
The  location  dependence  causes  the  “curse  of  dimensionality”  in 
stochastic  multiscale  simulations. 

□  Goal:  Consider  the  correlation  of  microstructures  between  different 
points.  Separate  random  variables  from  coordinates.  Dramatically 
reduce  the  dimensionality  of  the  stochastic  multiscale  input. 

Grain  structures 
of  a  nickel-base 
superalloy  turbine 
disk  having  dual- 
mi  crostructu  res 


(B.  Wen  and  N.  Zabaras,  2012) 
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Bi-orthognoal  Decomposition 


□  Start  from  realizations  of  the  microstructure  random  field  (A) 
varying  in  both  micro  (s)  and  macro-scale  (x):A(x,s,&>) :Jx5xQ->l 

□  Project  A  to  a  set  of  bi-orthogonal  bases 


A(x,s,«)  =  A(x,s)  + A(x,s,ro)  A(x,s):=f  A(x,s,  co)p(co)dco 

00 

=  A(x,  s)  +  £  4p,' 'P,  (s)O,.  (X,  co) 


i=l 

d 


A(X,  S)  +  X  \[PI'¥,  (S)®,  (X,  CO) 


i= 1 


□  The  inner  product  in  the  microstructure  domain  is  denoted  as 


□  The  inner  product  in  the  spatial  domain  is  defined  as 

{®, ,  ® j }  :=  Jx  (®,-  (x,  co)  ■  <D  (x,  co )}  dx 

□  The  orthogonality  conditions: 

^^J^^/entun^M/Vaj^nc^^^<amiaclakis^008^^<ouchmeshk^nc^^Zabaras^(M0^ 
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Bi-orthognoal  Decomposition 


□  By  minimizing  the  distance  between  the  Karhunen-Loeve  expansion 
and  the  random  field,  the  microscale  basis  can  be  computed  by 


%(*)  = 


l 


JPi 


A,0>, 


□  The  macroscale  basis  is  obtained  through  orthogonality  condition 

1  r  ~ 

o.(  x, &>)  =  —=!  A(x,s,&>  j'Tfsk/s 

V  P i 

□  The  eigenvalue  problem  in  K-L  expansion  can  be  defined  as 

p,>P,(s)  =  |vC(s,s')'P,(s')ds' 

□  The  covariance  matrix  is 

C(s,s')  =  |a(x,s,&>),  A(x,s’,&>)| 

n„,  n„ 


l  ",  "el  "int 

Or  in  discrete  form:  c(s,s')  =  — 


n. 


i=1  4=1 


if  realizations  of  the  random  field  are  given  as  A.(x,s,<».),  where  /  = 
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Second  Level  KLE 


□  The  random  field  A(x,s ,co)  is  decomposed  into  a  set  of  microscale 
modes  ^(s)  and  spatial-random  coupled  modes 

□  Q);  (x,&>)  is  still  high-dimensional.  By  assuming  independence  of 
these  macromodes,  the  random  variable  can  be  further  separated 
from  spatial  modes  via  a  second  level  KLE 

r- 

®,-  (x,  to)  «  0>  (x)+  ^^7 Yf  (x)<pf(ft>)  i=\,...,d 

M 

□  A,/  and  y/f  are  eigenvalues  and  eigenvectors,  respectively,  of  the 
covariance  matrix 

^  1  k= l 

□  The  dimensionality  of  the  reduced  representation  of  the  ith 
macromode  is  >)  .  The  total  dimensionality  r  of  the  reduced 
space  of  feature  A(x,s,<y)  is  then  the  sum  of  all  rt . 

d 

r  =  Y.r< 
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Polynomial  Chaos  Expansion 


□  After  performing  the  second  level  KLE  on  macro  modes  (x,&>), 
we  separate  the  randomness  from  spatial  coordinates.  The  random 
term  (p/  (a>)  can  be  mapped  to  well-shaped  distribution  by 
polynomial  chaos  expansion  (PCE): 


<p/H=2>  /**?(£>)) 


k 


The  choice  of  polynomial  T1;  of  random  variable  that  follows 
well-shaped  distribution  (e.g.  Gaussian  or  Uniform)  depends  on  the 
specific  distribution  of  C/  ■ 

□  In  the  current  work,  we  choose  C/  to  be  uniformly  distributed 
between  -1  and  1 .  Therefore,  Yf  are  Legendre  polynomials.  The 
coefficients  of  the  PCE  can  be  computed  by 
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Non-intrusive  Mapping 


□  (p j  ( co )  is  non-linearly  mapped  to  C/  through  cumulative  density 
function  (CDF)  and  based  on  the  given  realizations: 


»/(«>)  =  r/(f/),r/  =  F-,oF 


□  Two  assumptions  were  made  during  this  Bi-orthogonal  KLE-second 
level  KLE-PCE  process: 

>  Assumption  1 :  the  random  variable  co  is  independent  from  microscale 
and  macroscale  coordinates  (x  and  s). 

>  Assumption  2:  O., /  =  J  are  independent  from  each  other. 
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Procedure 


□  Training:  establish  the  mapping  between  microstructure  space  and 
reduced-order  space. 

□  Testing:  sample  in  the  low-D  space  and  reconstruct  real  microstructures 


•*  *  •  Jt  \1 

•  >•  A  l*V>  . 

■  \ 

*  «  4 

T 

1 

1 

1 

1 

Bi-orthoaonal  KLE 


i> 


_  a  _ 

A(x, s,  co)~  A(x, s)  +  Yj  (s)® , (x,  co) 

i= 1 

<  =» 


(a)  Multiscale  microstructure  sample 


Forging 

v  .(p )  +f  =o 

f-r0=ajubjp 

▼ 


Stochastic 

Multiscale 

Simulation 


(b)  Macro-random  coupled  modes 

Second-level  KLE  /V 

<D,.(x,m)  *<D,.(x) 

^7  +  Z  VT W  (*)<p/  («>) 


Iyf(x)r‘(,7) ; 

a 


(e)  Properties  distribution  (d)  Uniform  distribution 


(c)  Reduced  surrogate  space 
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Initial  Samples 


□  The  input  to  the  stochastic  simulation  is  a  set  of  ingot  samples  whose 
microstructures  at  different  points  are  random  but  correlated. 

□  The  input  ingots  are  resulted  from  pre-processes. 

□  To  obtain  the  random  input  ingots  having  correlated  microstructures, 
we  generate  them  through  a  random  deformation  process. 

>  A  set  of  workpieces  whose  surfaces  are  randomly  curved  are  pushed 
against  a  flat  die. 

>  The  microstructures  at  all  points  of  all  the  initial  ingots  are  assumed  to  be 
identical. 

>  The  resulted  microstructures  at  different  points  of  different  ingots  are 
distinct  due  to  the  random  shape  of  the  initial  ingots. 


X 
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Initial  Surface 


□  The  surface  of  initial  workpieces  are  described  by  a  degree  6  Bezier 
curve 


Rp(a,co)  =  0.5  1  +  YdPi{a>)9i{a) 


i= 0 


^0(a)  =  (l-a)6  (a)  =  15a4  (l-a)' 

cpx  (a)  =  6a (l- a)5  ^>5(a)  =  6a5(l  — a) 

(p2  (a)  =  15a2  (l-a)4  %(«)  =  «6 

(p2  (a)  =  20a3  (1- a)3  a-x/L 


0  0.2  0.4  0.6  0.8  1 


0  0.2  0.4  0.6  0.8  1 


□  p.  are  selected  to  be  uniformly  distributed  in  (-0.1, 0.1) . 

□  1000  initial  samples  are  generated.  The  resulted  microstructures  (or 
more  precisely,  texture)  are  the  input  to  the  stochastic  simulation. 
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Deterministic  Solver 


□  Each  Gauss  point  in  the  macroscale  is  linked  to  a  microstructure. 

□  The  deformation  of  the  microstructure  is  controlled  by  the  local 
deformation  gradient  in  the  macroscale  via  Taylor  hypothesis. 

□  The  mechanical  properties  of  the  point  is  evaluated  on  the 
microstructure  by  crystal  plasticity  constitutive  and  homogenized. 

□  Homogenization:  (t)  =  — f  TdF 


Mean 


Standard  deviation 


eajh  stress 

nti 


1000  initial  samples 
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Stochastic  Input 


□  Mesh  of  the  macroscale  workpiece:  10x6  quadrilateral  elements. 

□  Microstructure  representation:  grain  sizes  and  orientations  of  20 
grains. 

□  Uncertainty  source  A:  random  textures 


□  The  total  dimensionality  of 
the  stochastic  input: 

10x6x4x20x3  =  14400. 

□  The  first  3  macro  modes  are 
preserved  which  captures 
around  95%  of  the  total  energy. 

3 

A(x,s,&>) «  Z  \iPi' (s)°.  (x>  ®) 

i=\ 

□  The  expectation  of  macro 
modes  energy  is  identical  to 
the  corresponding  eigenvalue. 
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The  Second  Level  KLE 


□  The  macro  modes  are  further  decomposed  by  the  second  level  KLE. 

□  Each  decomposition  preserves  more  than  95%  of  the  total  energy  of 
the  macro  mode.  The  dimensionality  of  reduced  macro  modes  are 

rx  =2,  r2  =3,  r3  =3  ,  respectively.  The  dimensionality  of  the  reduced 
space  of  the  microstructure  feature  is  therefore  r  =  rx  +r2  +r3  =  8 . 
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PDFs  of  Reduced  Representations 


□  The  PDFs  of  reduced  microstructure  representations  q>/  obtained 
from  the  model  reduction  on  the  given  1000  samples  are  constructed 
by  kernel  density. 

□  The  PDFs  of  reconstructed  q>/  through  PCE  using  10000  random 
samples  from  the  uniform  distribution  U(-1 ,1)  are  also  plotted. 


Reconstruction  vs.  Test  Sample 


□  Take  an  arbitrary  microstructure  realization  and  project  it  to  the 
reduced  surrogate  space.  We  can  reconstruct  the  microstructure 
using  its  reduced  representation. 

□  The  first  step  is  to  reconstruct  its  macro  modes  through  PCE-KLE. 


texture 


Dim  ension 
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Reconstructed  Mean 
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Reconstructed  STD 
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Increase  Dimensionality  of  the  Reduced  Space 


□  Next,  we  increase  the  dimensionality  of  the  reduced  space  to  18,  so 
that  the  reduced  representations  capture  99%  of  the  total  energy  of 
each  macro  mode:  r  =  ri+r2+r3  =3  +  7  +  8  =  18 

□  The  reconstructed  macro  modes  of  a  single  texture  realization  is 
more  “realistic”  than  the  previous  8  dimensional  representation. 
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Convergence  with  Increasing  Dimensionality 


Increase  the  dimensionality  of  the  reduced  space  to  18  (3+7+8),  which  captures 
99%  of  the  total  energy  of  each  macro  mode. 


diff 


p  ~p 

1 18  dim  1  8  dim 


mean 
difference 


standard 

deviation 


std 

difference 


‘  18dim 
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Statistics  of  One  Point 


Equivalent  strain 
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Outline 


□  Introduction  and  motivation 

□  Uncertainty  quantification  at  a  single  material  point 

>  Investigating  mechanical  response  variability  of  single-phase 
polycrystalline  microstructures 

>  Investigating  variability  of  fatigue  indicator  parameters  of  two- 

phase  nickel-based  superalloy  microstructures 

□  Uncertainty  quantification  of  multiscale  deformation  process 

□  An  efficient  image-based  method  for  modeling  the  elasto-viscoplastic 
behavior  of  realistic  polycrystalline  microstructures 

□  Conclusion  and  future  research 


Cornell  University 

Materials  Process  Design  and  Control  Laboratory 


61 


Problem  Definition 


□ 


Goal:  Accurately  and  efficiently  investigate  effective  and  local  mechanical 
properties/response  of  polycrystalline  materials  based  on  realistic 
microstructure  image. 


□  Solution  strategy:  Green  function  method  in  combination  with  fast  Fourier 
transform  solving  governing  equations. 

□  Merit:  no  complex  meshing,  no  inversion  of  huge  matrix,  consider  both 
intergranular  and  intragranular  interactions,  take  image  as  input. 

□  New:  elasto-viscoplastic  constitutive,  application  to  fatigue  properties,  multigrid 
strategy. 

□  Local  mechanical  response  of  a  heterogeneous  medium  can  be  calculated  as 
a  convolution  integral  between  a  linear  homogeneous  reference  medium 
Green  function  and  the  actual  heterogeneous  field. 

□  Fast  Fourier  transform  is  introduced  to  reduce  convolution  integrals  in  real 
space  to  simple  products  in  Fourier  space. 
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Governing  Equations 


□  The  governing  equations  are  equilibrium  equations  and  boundary 
conditions.  An  separate  formulation  is  proposed  for  solving  crystal 
elasto-viscoplastic  problems,  where  elastic  and  plastic  responses  are 
computed  separately. 

Elastic  problem: 


□ 


□ 


Equilibrium  (stress  rate  divergence): 
Periodic  boundary  conditions. 
Elasto-plastic  problem: 

Equilibrium  (stress  divergence): 


IJJ 


=  0 


U,J 


=  0 


Periodic  boundary  conditions.  Incompressibility: 


4=0 


□  Solution  strategy: 

>  Solve  the  two  sets  equations  for  elastic  and  plastic  responses  separately. 

>  Represent  a  local  quantity  by  the  mean  plus  fluctuation. 

>  Reform  the  equilibrium  equation  using  Green  function  method. 

>  Transform  reformed  equations  to  Fourier  space  through  FFT. 

>  Update  strain  related  quantities  in  Fourier  space  then  transform  them 
back  to  real  space. 

>  Update  stress  related  quantities  in  real  space  following  constitutive  laws. 
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Elasticity:  Green  Function  Method 


□  Stress  rate  and  polarization  fields: 


a(\)  =  C‘, :  i:'  ( \ )  +  (p‘  (  \  ) 


Co  =  C‘ 


(pe  (x)  =  c(x)  -  C|:i : (x)  =  d(x)  -  Ce0 :  i:‘  (x) 


□  Equilibrium  equation:  £  .+0. .  =C‘v?r  +<pe  =0 

y,J  0 ykl  kl,j  t I, J  0 ijkl  k,lj  rij,j 

□  Solve  the  equilibrium  equations  using  Green  function  method: 

Assume  solutions:  v‘  =-|R3G^,(x-x’)/;(x')dx’  fem  (x)  =  (pemnn (x) 

Substitute  solutions  into  equilibrium  equations: 

-clm  L  GLji  (x -X (x’)dx'  +  JR3  Simf„ Z  (x  V(x  -  X  ')dx  ’ 1  =  0 

L  (x  -  x ')  +  4AX  -  x ')]  fn,  (x  ')dx '  =  0 


lf/me(x’)  is  arbitrary  -Ce0ijklGekmJJ(x-x')  +  SimS(x-x')  =  0 


As  a  result:  vf  =-|R3G;m(x-x>;.  .(x’)dx' 

LG:a(x-xv;v(x’)dx' 


\k = - 


^^^^bensohn^OC^ 
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Elasticity:  Green  Function  Method  and  FFT 


i 


Integrate  by  parts  and  assume  v‘  =  j\  G-mj{\ - x  V^.(x')dx' 
boundary  terms  vanish 

Va=|R3^(x-xV;.(x')dx’ 


^jCemGl^)  =  -Si 


□  Transform  to  Fourier  space: 

Governing  equation  in  Fourier  space: 

Solutions  in  Fourier  space: 

fft(v.) = % = i  FFm,) = t 


7  j  im  v  v/  t  mj 

Green  operator  in  Fourier  space:  fjH(£)  =  -£,£,<%(§)  <=>  Geikjl(%) 
Define  4  =  4,4 jCm  i=>  g;©  =  -(a-1  ) 


V^O  and  f4(0)  =  0 
Velocity  gradient  fluctuation 


|!(x)  =  l/2(v^.(x)  +  v;.(x)) 

^  (x)  =  1/2  ( v.j  (x)  -  Vjj  (x) ) 
Strain  rate  and  spin  fluctuation 
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Viscoplasticity:  Green  Function  and  FFT  Approach 


□  Stress  and  polarization  fields: 


g(x)  =  Cq  :  zp  (x)  +  (p^  (x) 


Cq  =(M 


P~  1 


(p^  (x)  =  o(x)  -  Cq  :  zp  (x)  =  g(x)  -Cp  :zp  (x) 

□  Equilibrium  equation:  6fy,y  -  ~0ijkrklJ  .  TijJ 

□  Equilibrium  equation  in  Green  function  form: 


=  CL,., £?., ,  +  tp? j  =  C^mvklj  +  cp? j  =  0 


~cmGLjj  (x  -  x')  +  4/(x  -  x’)  =  0 


□  Equilibrium  equation  in  Fourier  space: 

^jCpmjUGL^)=-sA  ^©=-(4)_1 


AP  —  Z  Z  Cp 
^ ik  ~  hlhj^Oijkl 


Solution:  FFT(vf)  =  vf  =i  FFT(v[k)  =  Kk  = 


v^(§)  =  r^(§)^(§) 
V^O  awJ  v/y  (0)  =  0 


^(x)  =  1/2  (v£(x)  +  v'(x)) 
<(x)  =  1/2(v5(x)-v;.(x)) 


Velocity  gradient  fluctuation  Strain  rate  and  spin  fluctuation 
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Elasto-viscoplastic  Solution 


□  The  total  velocity  gradient  combining  both  elastic  and  plastic  parts  at 
a  single  point  is:  Vv(x)  =  vv  +  Vve(x)  +  Vvp (x) 

□  Total  strain  rate:  stJ (x)  =  l/2(v  ;(x)  +  v;i(x)j 


□  Elastic  strain  rate  is  computed  following  a  Newton-Raphson  process 

r  dF 

Fi  («/  (x))  F(«c(x))  =  £'(x)  +  sp(x)  -  s(x)  =  0 


Ci(x)  =  ^(x)- 


d£‘ 


dF  dse  dsp  dzp  do  do  p 

+ - =  1,  + - =  L  +Mf  :d tlA  :C 


dze  dze  dte  H  do  do  dte 

□  Plastic  strain  rate.  gP(x)  =  g  (x)-£e(x)  cp(x)  =  £p(x)  — ^frace(£p(x)) 

□  Stress  and  stress  rate:  d(x)  =  Ce(x) :  £e(x)  o(x)  =  C^(o(x)) :  zp(x) 


□  Polarizations  for  the  next  iteration: 

<Pe  (x)  =  <i(x)  -Cg  :se  (x) 

q>p(x)  =  a(x)-Cp  \tp  (x) 
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Algorithm-Basic  Formulation 


□  Basic  formulation:  based  on  the  exact  expression  of  Green  function 
for  linear  elastic,  homogeneous  reference  material. 


For  each  time  step,  the  iterative  algorithm  can  be: 

(1)  At  the  beginning  of  the  1st  iteration,  give  an  initial  guess  to  the  total  velocity  gradient: 
JJW(x)  =  n Vv(x) .  Then  compute  initial  elastic  and  plastic  strain  rates,  J!£e(x)  and  J!^(x),  as 
well  as  stress  „+Xx)  and  stress  rate  Jfcfx),  using  elasto-plastic  constitutive  relations. 

(2)  Compute  the  elastic  and  plastic  polarization  fields,  V(x)  and  for  the  ith  iteration. 

(3)  Transform  the  polarization  fields  to  Fourier  space  via  FFT: 

y©=PFT(  V(x))  V©=PFT(  V(x)) 

(4)  Compute  velocity  gradient  fluctuations  in  the  Fourier  space  at  the  (i+1)th  iteration: 

:'$'©);  Vi;*0,  and“vr(0)=0 


1+ 1 
n+ 1 


V^©=^(F©:'^©);  Vi^O,  and“V^(0)=0 


(5)  Transform  velocity  gradient  fluctuations  back  to  the  real  space  through  inverse  FFT: 


i+ 1 
n+ 1 


W(x)=FFT'[  ™W©)  %W(*)  =  FFr'{  “W"©) 


(6)  Compute  strain  rate,  strain,  and  stress  fields  using  updated  velocity  gradient. 


(7)  Check  the  convergence  8  = 


i+ 1  „  i ,, 

e-  e 


i+ 1 


( 


div(i+1<r) 


1/2 


i+l 


m 


1/2 


If  not  converged,  repeat  steps  (2)  to  (7). 


i+l 


i+ 1 


o(0) 


Equilibrium 

error 


(H.  Moulinecand  P.  Suquet,  1998) 
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Microstructure  Representation  and  Discretization 


□  The  input  to  the  CPFFT  simulator  is  pixelized  image  with  orientation 
parameters  associated  with  each  pixel  (or  voxel  for  3D). 

□  Regular  grids  with  N1xN2  (for  2D)  pixels  or  N1xN2xN3  (for  3D) 
voxels  are  the  discretization  of  the  image  (bitmap  image). 

□  The  ith  coordinate  component  of  real  points  (pixel  or  voxel)  in  real 
space: 

L  L. :  length  of  the  microstructure  in  the  ith  direction 
N.\  Number  of  points  in  the  ith  direction 


x,  =  0,— ,2— 


N,  N, 

□  The  ith  coordinate  component  of  frequency  points  in  Fourier  space: 


<&"  = 


r 


v 


N, 


\ 


+  n 


—  ,  n  =  \,...,Ni  Ni:  even  ^  = 
L- 


N,+\  'll 


V 


+  n 

2 


n  =  l,...,Ni  Ni:  odd 
L; 


In  order  to  take  advantage  of  FFT,  the  point  number  in  each  direction  should 
be  taken  to  be  the  integer  power  of  2. 

Input  format 


PT1:  coord(l),  coord(2),  coord(3),  orient(l),  orient(2),  orient(3) 
PT2 :  coord(l),  coord (2),  coord(3),  orient(l),  orient(2),  orient(3) 
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Geometric  and  Crystallographic  Texture  Evolution 


□  The  material  grid  in  the  real  space  should  be  updated  after 

deformation  through:  .  . 

x,(Xp)  =  X,+(L.X,+v(X,))d/ 

□  The  computational  grid  after  deformation  should  stay  regular: 

xc  (Xc )  -  (l  +  Ld^ )  Xc 

□  For  convenience,  the  deformed  material  grid  is  usually  approximately 
assumed  to  be  coincident  with  the  regular  computation  grid  (single¬ 
grid  simplification): 

Xp  =Xc 

□  The  grain  orientation  is  rotated  according  to  the  total  spin  tensor: 

Ns 

»(x)  =  [n  +  cb(x)  -  ( bslip  (x)]  dt  ($slip  (x)  =  (x)/a  (x),  pa  (x)  =  antisym{  ) 
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Plane  Strain  Example 


□  Example:  A  FCC  aluminum  microstructure  containing  64  grains  (with 
random  orientation)  discretized  by  16x16x16  voxels  for  CPFFT  and 
16x16x16  elements  for  CPFEM. 

□  Random  initial  texture. 

□  Elastic  constants:  C11=110GPa,  C12=59GPa,  and  C44=26GPa. 


x  10 


□  Material  parameters:  r0=47,r1=86,6>0  =  550,6>=i6,/0=i 

□  Applied  velocity  gradient:  plane  strain: 


L  =  VV  =  10“3x 


0  0  0 

0  1  0 

0  0-1 


M'=Xo|i-a(X)0ma(X) 


f 


a= 1 


g“(X> 


m“(x) :  o(x) 


\ 


n— 1 


g“(X> 


Cornell  University 

Materials  Process  Design  and  Control  Laboratory 


71 


Micromechanical  Response 


□  Strain,  plastic  strain,  stress  fields,  and  texture  (single  grid): 

CVPFFT 


strain 


plastic 

strain 


stress 


texture 
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Homogenized  Mechanical  Response 


□  effective  stress  -  strain  curve  of  the  microstructure 


Convergence  with  Resolution 


□  An  example  of  the  same  microstructure  discretized  by  32x32x32 


voxels  is  also  conducted  to  show  the  convergence  with  finer  gird. 

EPFFT-16P  JL. 


250,- 


CEPFFT-16P 


CEPFFT-32P 


<  no  >  <  111  > 

CEPFFT-32P 


0.2  0.3  0.4 

Equiva lent  strain 


Plastic  strain 


- CEPFFT-16P 

- CEPFFT-32P 
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Multigrid  Strategy 


□  The  strain  related  fields  are  computed  through  FFT  on  computation 
grid. 

□  The  stress  related  fields  are  updated  according  constitutive  model  on 

material  grids.  \/(x„)  =  Z 


Real  velocity  grad 


o 


Constitutive  model 


Real  stress 


o 


Real  polarization 


Material  Grid 


Real  strain  rate 

V  J 

r  \ 

FFT  velocity  grad 

k _ / 

Real  velocity  grad 


O  invFFT 


Green  ^LJ^operator 


FFT  polarization 


O  FFT 


Real  polarization 


Computation  Grid 


^^^3hellec^^JVichel^HJVoulinec^anc^^uc|uet^0(^ 
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Mechanical  Response  and  Texture 


□  Effective  stress-strain  response  and  orientation  distribution  predicted 
by  CEPFFT. 
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A  Fatigue  Example 


□  Stress  -  strain  response  in  the  z-direction  under  cyclic  loading. 


The  same  INI 00  superalloy  as  used  in  the  FIPs  study. 


N. 


Mf  =2]m“(i)®m'(i) 


h 


Jh_ 

d: 


a  „  a 

a 

Tt  ~  Xx 

—  &A 

D", 
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nx- 1 


+  r2 


n ~ 


d: 
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Distribution  of  Grain  Level  FIPs 


□  The  distribution  of  FIPs  among  grains.  maxP  means  the  maximum 
FIP  of  an  individual  grain  and  aveP  means  the  average  value  of  the 
FIP  in  a  grain. 


(a)  CEPFFT,  16x16x16  voxels 

(b)  CEPFFT,  32x32x32  voxels 

(d)  CPFEM,  16x16x16  elements 

The  FIPs  predicted  by 
CEPFFT  are  more  spread  than 
CPFEM. 
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Grain  level  FlP/max  FIP 

(d) 


-maxPcyc 

-maxPFS 

-maxPmps 

-avePcyc 

-avePFS 

-avePmps 


-maxPcyc 

-maxPFS 

-maxPmps 

-avePcyc 

-avePFS 

-avePmps 
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Fields  of  FIPs 


□  Contour  plots  of  FIPs  over  the  microstructure. 


CEPFFT16P 


CEPFFT 32 P 


CPFEM  16E 


a-i«S 
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:  ::: 

0  CCSS 
OCOS 

ere  is 
ocoi 
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=^3 
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C  CC45 
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00015 
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c  cos 
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C  0045 
C  004 
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cccss 
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o  rei  5 


Pcyc  PFS 
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Computation  time  (s) 


Computation  Efficiency 


Plane  strain: 
Strain=0.1, 5<10-4 

100000 


Fatigue: 

1  complete  loop,  5<10-4 

100000 


1-Processor, 

CEPFFT, 

16xl6xl6-Pixel 

60-Processors, 

CEPFFT, 

16xl6xl6-Pixel 

1-Processor, 

CEPFFT, 

32x32x32Pixel 

240- Processor, 

CEPFFT, 

32x32x32Pixel 

1-Processor, 

CPFEM, 

16xl6xl6Element 

240-Processor, 

CPFEM, 

16xl6xl6Element 

164 

5 

1309 

17 

12687 

154 

554 

16 

4385 

57 

29026 

348 
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Outline 


□  Introduction  and  motivation 

□  Uncertainty  quantification  at  a  single  material  point 

>  Investigating  mechanical  response  variability  of  single-phase 
polycrystalline  microstructures 

>  Investigating  variability  of  fatigue  indicator  parameters  of  two- 
phase  nickel-based  superalloy  microstructures 

□  Uncertainty  quantification  of  multiscale  deformation  process 

□  An  efficient  image-based  method  for  modeling  the  elasto-viscoplastic 
behavior  of  realistic  polycrystalline  microstructures 

□  Conclusion  and  future  research 
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Conclusion 


□  Model  reduction  tools  are  developed  to  build  the  stochastic  input 
based  on  a  given  dataset  of  random  microstructures. 

□  Physics-based  deterministic  simulators  are  developed  to  evaluate 
material  properties/responses  according  to  underlying 
microstructures. 

□  Variability  of  material  properties/responses  induced  by  microstructure 
variation  is  studied.  The  uncertainty  quantification  is  performed  for 
microstructures  at  a  single  material  point  and  for  an  entire  workpiece, 
respectively. 
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Multiscale  Modeling  of  Superalloy  Systems 


Workpiece 


Grain  structure 


Dislocation  and  precipitates 


Macroscale  uncertainties 


Mesoscale  uncertainties 
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Topology 


Two-phase  features 


Orientation 


Model  Parameters:  CRSS,  etc. 


Location-specific  microstructures 


Process  parameters: 

•  Temperature 

•  Strain  rate 

•  etc. 


Initial  and  boundary  conditions 
(e.g.  Die  shape) 


Microscale  uncertainties 


Particle  size/volume  fraction 


Particle  shape 


Dislocation  confiouration 


Parameters:  APB  energy,  etc. 
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Discrete  Dislocation  Dynamics  in  Superalloys 

rm=200b,  APB  energy  density=200  mJ/m2,  volume  fraction  =10%. 


(a)  *  (b) 


Dislocations  in  the  nickel-base  supperalloy  containing  large  particles,  (a)  Simulated  dislocation 
configuration  interacting  with  precipitates  having  large  radius  (200b)  in  current  3D  DD 
framework,  (b)  Simulated  dislocation  configuration  in  2D  DD  framework  (Mohles,  2004). 


Future  focus:  massive  dislocation  dynamics  in  precipitation  hardened  alloys. 
Multiscale  linking. 
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Uncertainty  Quantification  with  Realistic  Polycrystalline  Microstructures 


□  Introduce  realistic  polycrystalline  microstructures  into  uncertainty 
quantification. 

□  Provide  more  reliable  prediction  to  material  properties. 

□  Account  for  higher  order  spatial  correlation  of  the  microstructure. 

□  Efficient  physics-based  deterministic  solver:  CEPFFT. 

□  Stochastic  microstructure  input: 

>  Pixelized  models 

•  Dimensionality  depends  on  resolution 

•  Need  preprocessing  on  input  data 

•  Sampling  of  new  samples 

>  Statistical  features 

•  Physical  constraints 

•  Microstructure  reconstruction 
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A  Preliminary  Example 


□  Grain  size  distribution:  Lognormal  distribution 

>  Mean  grain  size  (diameter)  =  0.1457 

>  Standard  deviation  =  0.0437 

>  Grain  number  =  60 

>  Assume  spherical/round  grains 

□  Random  crystallographic  texture  generation. 


32x32  pixels 


t=300s 

al,a2,a3  ~  U (-0.00 1,0.001) 


Deterministic  solver: 
CEPFFT 
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Mixture  of  PPCA 


□  The  mixture  of  PPCA  (MoPPCA)  is  a  combination  of  local  PPCAs  for  the 
data  that  is  approximately  piece-wise  linear. 

□  Two  steps  (the  two  steps  are  interactive): 

>  A  partition  of  data  into  classes. 

>  Fitting  local  PPCA  models  within  each  class. 

□  The  overall  model  distribution:  p(y)=Yt^Jp(y\j) 

□  Local  PPCA  model:  The  data  y  is  generated  by  a  linear  combination  of  latent 

variables  z:  y  =  Wz  +  ji  +  £ 

z~7V(0,I)  \i  \  mean  of  the  data  model  £~7V(0,<t2I) 

□  An  iterative  EM  algorithm  can  be  designed  to  find  the  optimization  of  all 
model  parameters: 

.,l  ,\V,rr; 

•  Generation  of  new  sample  y  requires  the  random  choice  of  a  local  model 
according  to  nj  ,  followed  by  sampling  as  in  single  PPCA. 

•  Two  ways  to  generate  new  data  samples  in  PPCA: 

>  Generate  from  the  marginal  distribution  of  data  y  directly: 

p( y)  =  N( y  |  Sampling  from  high-dimensional  data  space 

>  Generate  from  both  the  distribution  of  latent  variable  z  and  noise  £  : 

p( z)  =  N(z  1 0,1)  /?(£)  =  N(0, cr2MLl)  Sampling  from  low-dimensional  latent  space 
/  Determination  of  intrinsic  dimension:  Maximum  likelihood  estimation. 
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Variability  of  Mechanical  Responses 


Strain 

mean 


Strain 

STD 


Stress 

mean 


Stress 

STD 


Equivalent  stress-strain 
response  with  error  bars. 


>The  number  of  PPCA 
components  is  chosen  to  be  3. 
>The  ML  estimation  of  the 
dimension  of  the  latent  space 
is  3  . 


latent  space  data  space 
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Advanced  Methodology,  Property  Prediction  &  Material  Design 


□  Stochastic  Multiscale  Analysis  (Training): 


Parameters 

(Si.6i.Yi) 


Reduce  input 
complexity 


Model  reduction 
techniques 


Initial  Structures 


Compute  material 
response/properties 


Stochastic  multiscale  solver  (SPDEs) 


Material  Processing 


Crystal 

plasticity 


Dislocation 

dynamics 


Molecular 

dynamics/ 


Property  Variability 


□  SPDE  solver:  Bayesian  regression,  HDMR+ASGC. 

□  Prediction: 


A  new  set  of 
structural  and 
process  information 


Model  reduction 
and  classification 


PCA,  Manifold 
learning,  SVM,  x-mean 


Database 
obtained  through 
analysis  process 


□  Design  (Optimizing): 


Desired  properties 
and  initial  information 


classification 


PCA,  Manifold 
learning,  SVM,  x-mean 


Database 
obtained  through 
analysis  process 


Probabilistic  Prediction 


Bayesian  inference, 
Regression,  interpolation 


Design  and  feedback 


Bayesian  regression, 
^interpolation,  optimization. 


Distributions 
of  properties 
of  interest 


Proposed 
process  strategy, 
microstructure 
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Thank  you  very  much! 
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